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Abstract. A calculation of the pion balance functions in a thermal model is 
presented. The total result consists of resonance and non-resonance parts. A 
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1. Introduction 

In the last years the thermal (statistical) models of particle production in ultra- 
relativistic heavy-ion collisions turned out to be very successful in describing hadron 
yields [1-3], spectra [4], the elliptic-flow coefficients [5,6], and the HBT radii [5,7]. 
In view of this fact it is interesting to study further observables in the thermal 
approach and compare the results of the model calculations with the data. One 
particular example of such observables are balance functions [8,9], which describe 
correlation between opposite-charge particles in the rapidity space and are closely 
related to charge fluctuations [10,11]. 

The balance functions measured by the STAR Collaboration [12, 13] at RHIC 
are defined by the formula 



where (5) is the number of the opposite-charge pairs such that both members 

of the pair fall into the rapidity window Y and their relative rapidity is \y2 — yi | = S, 



.(6)) - (N ++ (6)} (N_ + (5)) - (N__(6)} 
(N + ) + (N.) 



}, (1) 
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Fig. 1. Contributions to the balance function from the neutral resonances (solid 
line) and from the non-resonance pions (dashed line), plotted as a function of 
the rapidity difference of the two pions. The temperature of the system is T = 
165 MeV, and the expansion parameters of the model correspond to the average 
transverse flow ((3±) =0.5. 



N+ is the number of positive particles in the interval Y, and other quantities are 
denned in an analogous way. 

The measurement [12] showed that the widths of the balance functions are 
smaller than expected from models discussed in Ref. [8] and significantly smaller 
than observed in elementary particle collisions. This problem was discussed by 
Bialas in Ref. [14], where the small widths of the balance functions were explained 
in the framework of the coalescence model [15]. In this paper we present the results 
of an alternative calculation of the 7r + 7r - balance function [16], which is based on 
a thermal model with resonances [4]. In a more recent paper by Pratt et al. [17], 
it was shown that the measured balance function may be well reproduced in the 
canonical blast wave model. 

In our approach, the 7r + 7r~ balance function has two contributions related to 
two different mechanisms of the creation of an opposite-charge pair. The first one 
(resonance contribution) is determined by the decays of neutral hadronic resonances, 
whereas the second one (non-resonance contribution) is related to other possible 
correlations among the charged particles. We assume that the second mechanism 
forces the two opposite-charge pions to be produced at the same space-time point 
with thermal velocities. On the other hand, the first contribution refers only to the 
decays of neutral resonances which have a 7r + 7r~ pair in the final state (we explicitly 
include K$, rj, vj p°, w, a, and /o). In this case the correlations among emitted 
pions are completely determined by the kinematics of the decays. The details of the 
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modeling of the non-resonance contribution and other technical remarks are given 
inRef. [16]. 
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Fig. 2. Balance functions for the pions in the thermal model calculated for four 
different centrality classes and compared to the experimental data of Ref. [12]. 
The normalization of the model curves was adjusted in each case and is listed 
near the plot labels. 

2. Results 

According to the discussion presented above, the 7r + 7r~ balance function can be 
constructed as a sum of the two terms 

B(S) = B R (5)+B NR (5). (2) 

The functions B R (S) and Bnr(S) resulting from our model calculation are presented 
separately in Fig. 1. The value of the temperature used in the calculation was T = 
165 MeV, and the expansion parameters of the model were fitted to the spectra of 
hadrons. This procedure yields the average transverse flow of 0.5c. One can observe 
that the widths of the two contributions are similar. The calculated total width, 
(S) = 0.66, turns out to be slightly larger than the experimental value for the most 
central collisions. The STAR result for the most central events (c = — 10%) is 
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(S) = 0.594 ± 0.019, for the mid-central (c = 10 - 40%) (S) = 0.622 ± 0.020, for the 
mid-periferal (c = 40-70%) (5) = 0.633±0.024, and for the periferal (c = 70-96%) 
(S) = 0.664 ± .029. Such dependence of the width of the balance function on 
centrality cannot be reproduced in our model by changes of the transverse flow 
within limits consistent with the single-particle spectra. 

In Fig. 2 our results are compared with the experimental values obtained for 
four different centrality classes. The normalization of the model curves was adjusted 
in each case, since we were not able to take into account, in a different way, the effect 
of a limited detector efficiency and acceptance. On the other hand, the kinematic 
cuts in pseudorapidity and transverse momentum were included exactly in [16]. The 
shapes of the model balance functions agree well with the data except for the most 
central case where the theoretical width is slightly larger. We note that dips of the 
experimental balance functions at very small values of 5 are caused by the HBT 
correlations. This type of effects is not included in our approach. We also note that 
the effects of the detector efficiency may influence the width of the balance function. 

As a conclusion we state that our simple thermal model with thermal and ex- 
pansion parameters fixed earlier by fitting the ratios of particle abundances and the 
transverse-momentum spectra gives a quite satisfactory description of the balance 
functions. 
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